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Abstract 

We investigate the resonant motion of neutral spin-^-fermions in a magnetic guide. A wealth of 
unitary and anti-unitary symmetries is revealed in particular giving rise to a two-fold degeneracy 
of the energy levels. To compute the energies and decay widths of a large number of resonances 
the complex scaling method is employed. We discuss the dependence of the lifetimes on the 
angular momentum of the resonance states. In this context the existence of so-called quasi-bound 
states is shown. In order to approximately calculate the resonance energies of such states a radial 
Schrodinger equation is derived which improves the well-known adiabatic approximation. The 
effects of an additionally applied homogeneous Ioffe field on the resonance energies and decay 
widths are also considered. The results are applied to the case of the 6 Li atom in the F = \ 
hyperfine ground state. 
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I. INTRODUCTION 



Magnetic traps are now being, used for many years in order to study and manipulate 
the properties of neutral atoms With the availability of efficient cooling techniques the 
controlled occupation of low lying quantum levels inside such traps became possible. This 
paved the way for the experimental exploration of several phenomena such as Bose-Einstein- 
condensation (see Ref. |2| and references therein) or the formation of degenerate Fermi gases 
P). Apart from studying such collective properties there is also a fundamental and vivid 
interest in the properties and the dynamics of single atoms. Its understanding is essential in 
order to identify systems for quantum information processing but it is also necessary to gain 
insights into trapping mechanisms in order to control the atomic motion. Theoretical studies 
of the quantum properties of neutral atoms in magnetic traps have been performed since 
the late eighties. Since then several field configurations have been in the focus of numerous 
prions, the q uachu P o, field Q, the whe trap BflQ or the noetic guide a nd the 
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Ioffe trap |g, ua, ■ However, only a few field configurations allow for stationary solutions. 
Bound states of spin- ^-particles trapped by a wire have been found analytically by Bliimel 

n n 

et al. |rjj and Vestergaard Hau et al. [6]. The latter authors employed a super-symmetric 
approach in order to derive the Rydberg series of bound states around an infinitely thin wire. 
Finite wire sizes were accounted for by introducing a quantum defect in the Rydberg formula. 
Burke et al. |7( were able to numerically derive similar series of energies for particles with spin 
up to 2 by using the finite element method combined with the application of multichannel 
quantum defect theory. However, in most cases atoms are trapped in meta-stable states 
which have a finite lifetime. The calculation of the corresponding resonance energies and 
lifetimes requires the application of sophisticated numerical methods. Bergeman et al. jj] 
have calculated dozens of resonances of spin-|-particles inside a three-dimensional magnetic 
quadrupole field by determining thephase-shift of scattered waves. A similar approach was 
pursued by Hinds and Eberlein li| which have studied the dynamics of neutral particles 
in a magnetic guide. Their results were reproduced by Potvliege and Zehnle who have 
performed a complex scaling calculation. In both publications the effect of an additionally 
applied homogeneous Ioffe field is also discussed but only a few resonances are studied. 
Although there are no principle obstacles to compute resonances there has always been 
put great effort in finding reasonable approximations which allow for an quick and easy 
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calculation of the resonance energies. Here the most prominent example is the so-called 
adiabatic approximation jj, where the spin is aligned along the local direction of the 
magnetic field. Depending on their electronic state the atoms are here divided into low- and 
high-field seekers which give rise to bound and unbound solutions. 

In this work we focus on the quantum dynamics of a neutral spin-|-particle in a magnetic 
guide. We also address the case of an additional homogeneous field being applied parallel to 
the guide (Ioffe field). Like Potvliege and Zehnle 1(| we employ a complex scaling approach. 
In contrast to their work we calculate hundreds of resonances for several values of the Ioffe 
field strength which enables us to analyze global features of the resonance spectrum. In 
order get further insights we develop an approximate description of so-called quasi-bound 
states. We derive a radial Schrodinger equation which is superior to the common adiabatic 
approach. In contrast to the latter our method becomes exact for high angular momenta 
and its quality does not depend on whether there is an additional Ioffe field or not. 

This work is organized as follows: In Section |n] we derive the Hamiltonian for a neu- 
tral particle possessing a magnetic moment /7 in a magnetic guide. We also account for 
an eventually applied homogeneous Ioffe field parallel to the guide. The resulting two- 
dimensional Hamiltonian depends on a single parameter involving both the gradient of the 
inhomogeneous field and the Ioffe field strength. It exhibits a plethora of unitary as well as 
anti-unitary symmetries. The effects of this intricate symmetry properties and in particular 
a resultant two-fold degeneracy of the resonance energies are discussed in Sec. IIHI In Sec. 
IIVI we briefly outline the numerical approach we pursue in order to obtain the resonance 
energies and decay width which is based on the complex scaling method. Sec. 0is devoted 
to a discussion of our results. We analyze the resonance spectrum for several values of the 
Ioffe field strength. Furthermore we investigate how the energies and lifetimes of the reso- 
nance states depend on their angular momentum. The properties of so-called quasi bound 
states are addressed. We provide a radial Schrodinger equation whose eigenenergies are very 
good approximations of the true resonance energies. A comparison to the commonly used 
adiabatic approximation is also performed. Finally the results are applied to the case of 6 Li 
in a magnetic field generated by a current carrying wire together with a homogeneous bias 
field. 
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II. THE HAMILTON! AN 



We consider the motion of a neutral spin-|-particle which carries the magnetic moment 
\x inside a magnetic quadrupole guide. The magnetic field is given by 



B{r) 



( bx ^ 
-by 

\ B ) 



(i) 



The parameters b and B determine the gradient and the field strength of an eventually 
applied homogeneous field along the z-axis (Ioffe field), respectively. Generally the Hamil- 
tonian for a neutral particle with the magnetic moment jl and mass M in a magnetic field 
reads 

1 



H 



3D 



2M 



{p 2 x +p 2 y +p 2 z )-fl-B 



(2) 



Since there is in our case (0) no explicit dependence on the ^-coordinate the momentum in 
^-direction p z is conserved. Employing plane waves for the wavefunction in ^-direction one 
obtains a two-dimensional Hamiltonian which describes the dynamics in the x — y plane. 

n 

Adopting atomic units |l7( and with /2 = — f 5 where S is the spin operator one yields 
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2D 



2M 



2 2 gM 
Px+Py + ^r- i hx(7 x -bya y + Ba z 



(3) 



Here g is the ^-factor of the particle. The components of S are given through the Pauli 



matrices <7j. Performing the scale transformation x 
omitting the bars one obtains 

1 



{^f-Yxi and p = 3 Pi and 



H 2D = H = - (pl+p 2 y + xa x - yo y + ^a z ) 



(4) 



with 7 = B (ifr) 3 - For 7 = the gradient b does not explicitly appear in equation (jljl. 
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Therefore we find here the energy level spacing to scale according to (-4r-) 3 - 



III. SYMMETRIES 



A. The case 7 = 



Analyzing the structure of the Hamiltonian (J2J) for 7 = we find 15 discrete symmetry 
operations. They are listed in table|lj Each symmetry is composed of a number of elementary 
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T PyI X yS\ 


TI xy S2 


TPxIxyS^ 


TP X Py I X y S2 



TABLE I: Symmetry operations of the Hamiltonian (j3J) for 7 = 0. Top part: unitary symmetries. 
Bottom part: anti-unitary symmetries. 

operations which are provided in table |HI All symmetry operations are either unitary or 
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TABLE II: Set of discrete operations out of which all symmetry operations of the Hamiltonian (J1J) 
can be composed. 

anti-unitary, the anti-unitary ones involving the conventional time reversal operator T. In 
spite of its simplicity the system (J3J) possesses a wealth of symmetry properties. The algebra 
of the underlying symmetry group possesses a complicated structure some features of which 
will be discussed in the following. The operators E 2 , E y and E 2 generate a sub-group obeying 
the algebra [Ej, E^] = 2i eijk^k reminiscent of angular momentum operators. In addition one 
finds Ef = 1. Interestingly these quantities act on both real and spin space. 

Apart from the discrete symmetries there is a continuous symmetry namely A 2 = L z — S z 
which was already pointed out by Eberlein and Hinds in Ref. jjj. The operator A z obeys 
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the eigenvalue equation 



A z \m) = m \m) (5) 

with the half- integer quantum numbers m. For spin-|-particles the states \m) read in the 
spatial representation 

/ a T gi(m+§)A 
I a l e *(m- 3 )0 J 

The symmetries listed in table H] together with A 2 allow several sets of commuting operators. 
For our investigation we choose the set consisting of H, E 2 and A z . Calculating the eigenvalue 
of a A 2 eigenstate with respect to the S z operator leads to 

(Q,T e i(m+i)(<H7r) \ 
= e i(m+ 2 ) "|m) = /dm). (7) 

Therefore in case of m + | being odd we find k = —1 whereas we obtain k = 1 for even 
values of m + |. 

Exploiting the discrete symmetries one can show that there is a two-fold degeneracy of 
any energy level. This degeneracy is revealed as follows. The operations S 2 and P y a x obey 
{S 2 ,P 3/ cr a; } = 0. Let \E, k) be an energy eigenstate and at the same time an eigenstate of 
S z with 

E 2 \E,k) = k\E,k) (8) 
and k — ±1. Employing the above anti-commutator one obtains 

T^ z P y a x \E, k) = -P y ajl z \E, k) = -nP y a x \E, k) (9) 



The state P y a x \E, k) can be identified with \E, —k). Hence, as long as n ^ Vt^ there is 
always an orthogonal pair of states possessing the same energy namely \E, k) and \E, —k). 
We have to emphasize that there occur no further degeneracies in the system. In principle 
one could think of performing the above calculation repeatedly but now substituting P y a x 
with any operator listed in table H] which anti-commutes with E z . It turns out that all of the 
resulting states generated by this scheme are either superpositions of \E, k) and \E, —k) or 
differ only by a phase factor from one of these states. 
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P T <?* 

1 x ± xy u 2 




TPy 


TP x a z 


TP X Py I X y S2 


TI X yS2 



TABLE III: Discrete symmetry operations of the Hamiltonian @ for a finite Ioffe field strength 
B, i.e. 7 ^ 0. Top part: unitary symmetries. Bottom part: anti-unitary symmetries, symmetries. 

Since S 2 and A z commute, i.e. [E z , A J = 0, one can construct eigenstates with respect to 
both of these operators: \E, m, k). Letting P y a x act on these states together with considering 
the anti-commutator {P y a x ,A z } = one finds 

A z P y cr x \E, m, k) = —P y a x A z \E, m, k) = —m P y a x \E, m, k) . (10) 

Hence we can identify P y a x \E, m, k) with \E,—m,K). Therefore the pairs of degenerate 
states do not only have opposite k quantum number but also opposite m- values. A detailed 
discussion of degeneracies of spin-i-particles can be found in [12|. 

B. The case 7 / 

For the conserved quantity A z persists but only 7 of the discrete symmetries for 

7 = remain. They are listed in table ILTT1 These operations form a non-Abelian algebra. In 
contrast to the group operations listed in table |J there are no two anti-commuting operators. 
Hence it is not possible to construct pairs of degenerate energy eigenstates as discussed above. 
The operation P y a x which we have used to construct the degenerate pairs of states for 7 = 
now has the following property 

P v a x H(j) P y a x = H(-i), (11) 

i.e. the sign of 7 in the Hamiltonian (jlj) is changed. Hence, the application of P y a x reverses 
the sign of the homogeneous Ioffe field. We emphasize that even for 7 7^ the operators H, 
H z and A z form a commuting set of operators. 

IV. NUMERICAL TREATMENT 

The Hamiltonian (jlj) does not support bound states . In order to calculate the energies 
and decay widths of the scattering wavefunctions we employ the complex scaling method 
together with the linear variational principle. 
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We perform a rotation of the spatial coordinates into the complex plane: Xi — *> Xit %e 
which yields the 'rotated' Hamiltonian 

H(0) = 1 - (pie-™ + ple~™ + xe i6 o x - ye*a y ) . (12) 

Unlike the bound states only continuum states are affected by the complex scaling procedure. 
The scattering states are rotated into the lower half of the complex plane and in particular 
resonances correspond now to square integrable functions. Since for 9 ^ the complex 
scaled Hamiltonian ()12j) is not Hermitian one encounters complex energies: 

e = E-i T -. (13) 

Here E is the energy and V the decay width of the resonance (the lifetime r of a quantum 
states is connected to its decay width by r = T _1 ). In general the energy E will depend on 
the value of the rotation angle 9. Continuum states are rotated by 29 into the lower half 
of the complex plane whereas energies of resonance states - once revealed - are independent 
of the angle 9. A detailed discussion of the complex scaling procedure can be found in Ref. 
Q| and Refs. therein. Since the complex scaling procedure transforms resonance states 
into square integrable wavefunctions one can apply the linear variational principle in order 
to compute the wavefunctions and resonance energies. We utilize an orthonormal basis set 
of the form 

\m, n, s) = \m, n) \s) (14) 
where the functions \n, m) are products of harmonic oscillator eigenfunctions 

\m,n) = 7 J§2L ] , e-K^W)ff m Hn (^) (15) 

V2 m+ri 7rm! n\ 

The parameters uj x and u> y are non-linear variational parameters which can be adapted in 
order to gain an optimal convergence behavior of the numerically obtained wavefunctions. 
To cover the spin space dynamics we utilize the spinor-orbitals \s) = ||) and |f), respectively. 
Performing the linear variational principle the solutions \E) of the stationary Schrodinger 
equation are expanded in a finite set of the functions (fTHjl : 

\E) = J^c mn \m,n,s) (16) 
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From our knowledge of the symmetry properties we can choose a more specific appearance 
of this expansion. Requiring the states \E) to be eigenstates with respect to £ 2 we make 
use of the properties 



E z \m,n, t) = (-l) m+n \m,n, T) 
E z \m,n, i) = (-l) m+n+1 \m,n, I) . 

Thus the expansions for the two k subspaces read 

|£,« = +1) = ^2 a mn\m,n, t) + b mn\m,n, I) 



m+n=even 



(17) 
(18) 

(19) 

_ (20) 

m+n=odd m+n=even 

Since [A 2 , E 2 ] = one could also demand the states to be a priori eigenfunctions of A 2 as 
well. However, our basis set is not particularly well suited to construct eigenfunctions to 
A 2 . We therefore abstain from putting this constraint on the above expansion. Neverthe- 
less, the resulting diagonalization of the Hamiltonian matrix (see below) guarantees that 
our numerically obtained energy and E 2 -eigenfunctions are also eigenfunctions to A 2 . The 



m+n=odd 

\E,k=—1) = ^ a m n\m,n,'\) + ^ b mn\m,n,l) 
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FIG. 1: Dependence of the resonance energies (upper part) and decay widths (lower part) on the 
complex scaling angle 9 as observed in the numerical calulations. The stationary region is located 
in between the dashed lines. 

determination of the optimal expansion coefficients a mn and b mn gives rise to a complex 
algebraic eigenvalue problem H (8) v(9) = e (8) v (8) , where H is the matrix representation 
of the Hamiltonian (|12)1 . The vector v contains the coefficients a mn and b mn . The eigenval- 
ues of H (9) are computed for < 8 < j. Due to the finite size of the employed basis set 



resonances become not independent of the complex scaling angle 6 but exhibit a stationary 
behavior within a certain ^-interval as can be seen in figure H 



V. RESULTS 



In this section we discuss the results obtained via the above-described complex scaling 
approach. We present the resonance energy spectrum and study the dependence of the decay 
widths on the eigenvalue of A 2 . Additionally we discuss how the lifetimes are affected if a 
homogeneous Ioffe field ist applied. Furthermore a radial Schrodinger equation is derived 
that describes quasi-bound states in the magnetic guide. All results are given for b = M = 
| = 1. In order to obtain the resonance energies and decay widths for different parameters 

2 

one has to multiply the latter by the factor jj f-jn 3 . The section concludes with a discussion 
of how the results apply to the experimentally important case of a trapped 6 Li atom in the 
F — | hyperfine ground state. 



A. Energies and decay widths of resonance states 

In figure 121 we present the energies and decay widths of the resonances of a spin-|-particle 
in the magnetic guide. For 7 = one observes a regular pattern for the arrangement of the 
complex energies, i.e. they are located along lines in the semi-logarithmic representation. 
This implies that the decay width essentially decreases exponentially with increasing energy 
along the mentioned lines. The decay width of the lowest (and here also broadest) resonance 
is r^Y° = 0.3738. Increasing 7 results in a distortion of the regular pattern. One observes 
the formation of vertical lines together with a regrouping of the resonances into pairs for 
larger 7 (see in particular figure 121 for 7 = 4). Here states with large decay widths are 
affected more strongly than long lived states. Overall we encounter a decrease of the decay 
width if 7 increases. Thus as it is well-known an additional Ioffe field has a stabilizing 
effect. For the decay widths of the lowest resonance state our calculation yields the series 
r£rj = 0.0648, r^= n 2 = 3.08 • 10" 3 , r^= n 4 = 3.38 • 10~ 7 and r m = n 5 = 1.04 • 10~ 9 . These values 
suggest an exponential decrease of the decay widths if 7 increases. This agrees with the 
results of Sukumar and Brink who have estimated an exponential increase of lifetimes with 



growing Ioffe field strength by an analytical calculation 
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FIG. 2: Energies E and decay widths T (logarithmic scale) for resonances in a magnetic guide 
for k = 1. An increasing Ioffe field strength leads to a global decrease of the decay widths and 
therefore has a stabilizing effect. At the same time the resonance pattern becomes distorted and a 
regrouping of the resonances into pairs is observed. 



In Sec. IIHI we have pointed out the degeneracy of the k — 1- and k = — 1-subspaces. 
For 7 ^ these degeneracies are expected to split up since the symmetry properties of the 
system changes. Figure H3 shows the resonances for 7 = 1. States belonging to the k = 1 or 
K — — 1 subspaces are indicated by a dot or cross, respectively. Each state is accompanied by 
its degenerate partner for 7 = 0. The respective energy splitting decreases with decreasing 
decay width. In the following we will find that a small decay width is associated to the A z 
eigenvalue property of a resonance state: States with a small value for T are localized far 
away from the center of the guide. Here the quadrupole field dominates the homogeneous 
Ioffe field and determines the appearence of the spectrum. Therefore these states become 
less affected by the Ioffe field which results in a decrease of the energy splitting. 

Let us compare our results to those given by Hinds and Eberlein in Ref. [21]. First of all 
we remark that we have studied a significant larger number of resonances. In table IIVI we 
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FIG. 3: Resonance energies and decay widths in the k = 1- or n = — 1-subspaces at 7 = 1 (linear 
scale plot). Due to the finite Ioffe field formally degenerate pairs split up. States with negative 
quantum number m are shifted towards higher resonance energies and vice versa. The splitting 
decreases with decreasing decay width. 





m= \ 


m = | 




e r 


e r 





1.2937 0.3738 


1.7591 0.1122 




1.32 0.34 


1.765 0.11 


1 


2.1404 0.3788 


2.5144 0.1444 




2.125 0.34 


2.52 0.15 
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2.8438 0.3726 


3.1711 0.1594 




2.81 0.34 


3.17 0.18 



TABLE IV: Comparison between our results (bold face) and those obtained by Hinds and Eberlein 
for the first 3 energetically lowest resonances in the m = ^ and m = | subspaces. 

have listed the resonance energies and decay widths of the first 3 states inside the m= \ and 
m = I subspaces. Both results agree within a few percent. The difference might originate 
from the^complicated procedure employed by Hinds and Eberlein to locate the resonance 
energies 
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B. Dependence of the resonance energies on the eigenvalue of A z 



In this section we analyze how the energies and decay widths of the resonances are related 
to the eigenvalues of the angular momentum operator A z . Although A z is a conserved 
quantity we have employed a set of basis functions which do not respect this fact. Thus 
we need to calculate the matrix element m = (E, k\ A z \E, k) to receive the A 2 -eigenvalues 
of the resonances. In order to do this one has to respect the non-Hermitian character of 
the Hamiltonian (|T2j) which requires a complex symmetric scalar product. For a detailed 
discussion see Ref. Figure |U depicts the energies E and the A z eigenvalue of the 
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FIG. 4: Resonance energies and A 2 eigenvalues of resonances for 7 = (left picture) and 7 = 1 
(right picture). For 7 = the energies of the resonances form a symmetric pyramid-like pattern. 
For finite Ioffe field strength the pattern becomes asymmetric. 

resonances. For 7 = one observes each state to possess a degenerate counterpart, i.e. a 
state with opposite m. This results in the formation of a symmetric pyramid-like distribution 
where the maximum A z eigenvalue depends approximately linear on the resonance energy E. 
With increasing 7 the resonance energies in general are shifted towards larger values (this 
can hardly be seen in figure 0] for 7 = 1). Thereby, states with positive m acquire a larger 
energetical shift than states with negative values of the quantum number m. Due to this 
asymmetric energy shift the distribution becomes asymmetric (with respect to m — > — m), 
too. Here energies of states having the same k quantum number form continuous, nearly 
horizontal, lines whereas we observe the resonance energies for 7 = to be arranged on 
broken horizontal lines. In subsection IV Dl we will show that in the limit of ^ — > one 
finds a pattern of equidistant straight horizontal lines formed by the resonance energies. 
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Subsequent lines belong to states with opposite S z quantum number. With increasing m 
values the states become less affected by the Ioffe field. Their wavefunctions become located 
farther away from the center of the guide. In this region the strength of the linearly increasing 
quadrupole field outweighs the effects of the Ioffe field. 
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FIG. 5: Decay widths and A z eigenvalues of resonances for 7 = (left picture) and 7 = 1 (right 
picture). The decay width decreases exponentially with increasing modulus of m. With increasing 
Ioffe field strength one observes an overall decrease of the decay widths and a widening of the 
distribution. 



In figure the dependence of the decay widths on the A z eigenvalue is presented. Here 
we encounter an exponential decrease of the decay widths with increasing modulus of m. 
Overall the decay widths decrease if 7 increases and at the same time the distribution 
becomes broader. 



C. Quasi-bound states for 7 = 



Although a spin in our field configuration does not possess bound states there exist 
resonances with very long lifetimes, i.e. very small decay widths. The discussion of these 
quasi-bound states is the subject of the present section. 

Transforming the Hamiltonian (j3J) to polar coordinates (p, <p) and exploiting the conser- 
vation of A z yields 

,2 



H m = (m\ H \m) 



dp 1 
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pdp 
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+ P&x 



(21) 
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Performing the spin-space transformation 
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and introducing the spinor = p 2 |$) the corresponding Schrodinger equation becomes 
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(23) 



The transformation (|22j) diagonalizes the pB interaction term but leads to off-diagonal 
elements in the angular momentum term. In the limit of large m this coupling between the 
up and down components of |$) can be neglected: 



m + ma x 
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which yields 



d 2 
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(24) 



(25) 



By construction this radial Schrodinger equation (|23j) does not couple the up- and down- 
component of the spinor wavefunction |$). The lower component is unbound since the 
corresponding effective potential is V~ (p) = ^4 — \p- 111 contrast to this the potential for 
the upper component is V + (p) = |4 + \p and therefore bound solutions are allowed. We 
will refer to the corresponding states as the quasi-bound states \x) = (x (p) j 0) which are 
described by the Schrodinger equation 



m 



-X + —fX + PX 
P 



EqbX- 



(26) 



Here a prime denotes the derivative with respect to p. In order to solve equation (|2fijl we 
have utilized the FEMLAB software package wich employs the finite element method for 
solving differential equations (^(J. 

In table |V] the resonance energies E obtained from the complex scaling calculation are 
compared to the approximate energies E q b resulting from equation (|25J) . Even for m = \ 
we find a remarkable good agreement between E and E q \, although the validity of equation 
(f23j) is not justified since the off-diagonal coupling terms are of the same magnitude than 
the diagonal ones. With increasing m the discrepancy of E and E qb decreases as suggested 
by (|2*4*|) . The surprisingly good quality of the approximation can be explained by looking at 
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1 


2 


3 


4 


5 




E 


1.2937 


2.1404 


2.8438 


3.4688 


4.0419 


4.5769 


m= \ 


E q b 


1.2829 


2.1287 


2.8324 


3.4579 


4.0315 


4.5669 




% 


0.83 


0.55 


0.40 


0.31 


0.26 


0.22 




E 


3.2660 


3.8412 


4.3797 


4.8889 


5.3741 


5.8392 




Eqb 


3.2542 


3.8256 


4.3612 


4.8684 


5.3522 


5.8162 




% 


0.36 


0.41 


0.42 


0.42 


0.41 


0.39 




E 


4.7847 


5.2629 


5.7230 


6.1673 


6.5977 


7.0156 


m = f 


Eqb 


4.7805 


5.2577 


5.7169 


6.1602 


6.5897 


7.0068 




% 


0.09 


0.10 


0.11 


0.12 


0.12 


0.13 




E 


6.0956 


6.5208 


6.9344 


7.3375 


7.7311 


8.1160 


m = ^ 


Eqb 


6.0934 


6.5181 


6.9313 


7.3341 


7.7273 


8.1118 




% 


0.04 


0.04 


0.04 


0.05 


0.05 


0.05 




E 


7.2790 


7.6689 


8.0505 


8.4245 








Eqb 


7.2774 


7.6671 


8.0486 


8.4224 


8.7891 


9.1492 




% 


0.02 


0.02 


0.02 


0.02 







TABLE V: Comparison of the resonance energies to the approximate energies E q b obtained from 
equation (|25j) . The first six resonance energies E for 5 selected values of the quantum number m 
are provided. The rows labelled by '%' show the relative difference between E and E q b in percent. 

how the bound solution \x) = (x(p) > 0) T and the unbound wavefunction \() = (0,C(p)) T 
are coupled by the Schrodinger equation PHj) : 

?<xl^flC> = ? [dpx*(p)- 2 C(p) (27) 
2 p l 2 J p z 

From this expression one recognizes that transitions from the bound state \x) to the unbound 
state |C) are going to happen essentially at the center of the guide. However, since m can 
only adopt half-integer values the centrifugal barrier in equation (J25|) always persists. Thus 
both wavefunctions \x) and |C) vanish for p — > 0. However, this is the only region where the 
operator contributes significantly. Hence we have (x| % \() ^ and thus there is only 
a small coupling between the bound and unbound solution. Therefore the resonance states 
are very well described by the solutions \x) of equation This also explains the long 
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lifetime of states with high m quantum numbers. Here the large angular momentum barrier 
prevents the coupling between the bound and unbound channel. The particles are then 
located between the classical turning points of the potential Vqb(p) = + \p at a distance 
of approximately p ni \ n 



Performing a harmonic approximation of the potential 



10 

LLI 5 





m = 20.5 



m = 2.5 
m = 1.5 
m = 0.5 



FIG. 6: Energies E\ r of the lowest resonances of each m subspace calculated by using equation 
()28|) . Their accuracy is higher than the resolution of the plot. 



V + (p) around its minimum at p m i n yields a useful expression for the energy of the lowest 
resonance in each subspace 



(28) 



This equation provides a very good approximation to the energies. The calculated resonance 
energies differ from the energies E\ r by at most 0.05%. Figure El shows the energy spectrum 
obtained from equation (f2*%|) . For large m values the level spacing scales according to m~3. 



D. Quasi-bound states for 7 7^ and adiabatic approximation 



In the previous section we have shown the existence of quasi-bound states in case of a 
vanishing Ioffe field. We will now show that similar states exist also if a Ioffe field is present. 
We start our investigation with the Hamiltonian 

,2 



H m = (m\ H \m) = ~ 







Id (m+ \a.y 

■) 2 " ~7T + + P°* 

op z p op p z 



(29) 
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Introducing the spinor wavefunction = p 2 |$) and applying the unitary transformation 



a 




(30) 



with a 



the corresponding Schrodinger equation becomes 



d 2 m 2 ma 7 a 
H s- H ^-cr^ + -<7 2 + 



<9p : 



a 



4i 2 



— 1 



a 



7 <9p ^7 3 



ft, + 



ma 



IP 



|$) = (31) 



Like without Ioffe field the unitarian (J3*U|) transforms the /2i?-coupling term into a diagonal 
form. However, unlike the transformation (J22j) U in equation (|3Lip depends explicitly on the 
coordinate p. Therefore the transformation of the derivative results in additional terms: 



d 2 

U+ W u 



d d 2 



U + U" + 2U + U'— + 



dp dp 2 ' 

Neglecting the off-diagonal coupling terms in equation (|3*Tf yields: 



(32) 



d 2 m 2 ma 7 a 4 

+ -7T + -5-<T a + -<X Z + 



9p ; 



/0- 



fr 



a 



472 



1$) = E gb |$) 



(33) 



Here similar to equation pHjl only the upper component of the spinor | <£>) is bound. It obeys 
the Schrodinger equation 



O 2 



m /-= 

+ —r + Vl +P 2 + 



dp 2 



+ 



r 



X = E qb x- 



(34) 



Table IVII compares the resonance energies E obtained for 7 = 5 to the eigenvalues E q b 
calculated by solving the scalar radial Schrodinger equation ([34)1 . Apart from m = —0.5 
the agreement is excellent with discrepancies lower than 0.05 %. This is again the result of 
a localization of the particle wavefunction away from center of the guide which is the only 
region where the off-diagonal coupling terms are remarkable. However, for m = —0.5 and 
7 = 5 the effective potential of the Schrodinger equation (jMj) does not possess a centrifugal 
barrier. Here the wavefunction has non-zero contributions in the vicinity of the center 
of the guide. Hence the off-diagonal coupling terms of equation (J31|) become important. 
Considering the fact that the equation (J34j) here becomes certainly invalid the energies E 
and E qb still agree surprisingly well (see table IVljl . 

We now compare the results of the approximate Schrodinger equation (jMj) to those one 
would obtain within the so called adiabatic approximation. In this picture one assumes the 
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1 


2 


3 


4 


5 




E 


2.8163 


3.3978 


3.9412 


4.4547 


4.9437 


5.4123 


m = -\ 


E q b 


2.8471 


3.4272 


3.9691 


4.4813 


4.9692 


5.4369 




% 


1.08 


0.86 


0.70 


0.59 


0.51 


0.45 


I = 


E a d 


2.8387 


3.4205 


3.9635 


4.4764 


4.9649 


5.4329 




% 


0.79 


0.66 


0.56 


0.48 


0.43 


0.38 




E 


4.4663 


4.953 


5.4199 


5.8696 


6.3045 


6.7262 


m = ± 


Eqb 


4.4648 


4.9514 


5.4180 


5.8677 


6.3024 


6.7240 




% 


0.03 


0.03 


0.04 


0.03 


0.03 


0.03 


I = 6 


E a d 


4.4957 


4.9797 


5.4442 


5.8921 


6.3253 


6.7457 




% 


0.65 


0.54 


0.45 


0.38 


0.33 


0.29 




E 


6.8534 


7.2540 


7.6453 


8.0280 


8.4028 


8.7703 


m = f 


Eqb 


6.8523 


7.2528 


7.6440 


8.0267 


8.4014 


8.7688 




% 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 


I = 16 


E a d 


6.9044 


7.3027 


7.6918 


8.0725 


8.4455 


8.8114 




% 


0.74 


0.67 


0.60 


0.55 


0.51 


0.47 



TABLE VI: Comparison of the resonance energies to the approximate energies E q b and adiabatic 
energies E a d for 7 = 5. The first six resonance energies E for 5 selected values of the m quantum 
number are given. The rows labelled by '%' show the relative difference between E and E q b or 
Ead, respectively, in percent. 



projection of the atomic spin onto the local direction of the magnetic field to be conserved. 



B 



with ms being 



Thus the coupling of the magnetic moment to the field reduces to gfiB^s 
the projection of the spin onto the local field direction. In case of a spin-|-particle in the 
magnetic guide the corresponding Hamiltonian becomes 

1 



H. 



ad 



Px+Py±Vi+ x +y z 



(35) 



having employed scaled coordinates (see Sec. HIj) . Considering only the positive sign (which 
allows bound solutions) and introducing the wavefunction \^ a d) = p ^ \®ad) the correspond- 
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ing Schrodinger equation becomes 



d 2 i 2 
+ 



dp 



+ Vi 2 Tf 



\<&ad) = E ad \$ ad ) 



(36) 



Here I is the quantum number of the operator L z which is conserved due to the rotational 
invariance of the system around the z-axis. Note that unlike m the quantum number / is 
integer- valued. Table IV II shows a comparison of the adiabatic eigenvalues E ad to the exact 
resonance energies as well as the energies of the quasi-bound states E q b. The quantum num- 
bers / and m are chosen such that I = 771+5. Only for I = or m = — \ the adiabatic 
energies E a d are in better agreement to the exact ones than the quasi-bound energies E q b. 
One has to note that in contrast to E ad the energies E q b become exact in the limit of high m 
quantum numbers. Moreover we have to emphasize that only the Schrodinger equation p4|) 
yielding the quasi-bound states reproduces the correct degeneracies of the system. In con- 
trast to this the corresponding adiabatic equation shows a two-fold degeneracy of the states 
\l) and |— I) for any value of 7. We now investigate how the different approximations perform 




FIG. 7: Comparison of the quasi-bound energies E q b and adiabatic energies E a d to the exact 
energies E. The figure shows the energies of the lowest four states in the m = =^ and I = 8 subspace, 
respectively. The discrepancy between E and E q b is hardly visible throughout the complete 7- 
interval shown. 



for different 7 values. In figure (J7J) the energies obtained from each of the three methods 
(exact, quasi-bound, adiabatic) are depicted for the 4 energetically lowest resonances in the 
m = W and I = 8 subspace, respectively. One observes a remarkable agreement of E and 
E q b throughout the complete 7-interval. In contrast to that severe discrepancies between the 
adiabatic and exact energies are revealed for small values of 7. This shows the extremely 



20 



good performance of the quasi-bound approximation independently of the value of 7. This 
makes equation (J53J) the correct choice to calculate approximate eigenvalues. 

Figure Q also shows that if 7 becomes large the adiabatic approximation begins to 

2 

perform well. In the limit of ^ < 1 equation (J36j) can be further simplified. A series 

2 

expansion of the potential term up to first order in ^ yields 



If 



2 _ 1 

4 



dp 2 



P 1 



P j_ 

2 7 



\®ad) = Kd \$ad) ■ 



(37) 



This is the Schrodinger equation of a radial harmonic oscillator. The corresponding eigenen- 
ergies are given through 



E 



7 



ni 



+ 




2 7 



(2n+ \l\ + 1) 



(38) 



Introducing the frequency oj ■ 
the formula 

E„ 



^ together with the substitution / = m + | one arrives at 



7 



u I 2n 



(39) 



Figure |H1 shows a plot of the energies E nm . The corresponding S 2 -eigenvalues are calculated 
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FIG. 8: Graphical representation of the eigenenergies (|39|) . The values of the quantum number 
k are calculated using equation (J7J). A pattern of equidistant horizontal lines with alternating k 
values is formed. 



using equation (jJJ). The resultant pattern is similar to the one we have already observed in 
figure 0Jd. We find alternating equidistant horizontal lines of energies belonging to states of 
the two different S^-subspaces. 
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Resonance energies and decay widths of 6 Li in the 2 2 S 1 , F = ^ hyperfine ground 

2 Z 



In order to experimentally prepare a magnetic guide one can superimpose the magnetic 
field of a current carrying wire by a so-called bias field oriented perpendicular to the current 
flow The resulting magnetic field possesses a line of zero field strength at a distance 
Po = 2k°b P arane l to the wire. A series expansion of the field around po yields 



B 



Bb 
Po 



( „ \ ( _ x 2 + 2xy + y 2 \ / „f n *_1^\ \ 



X 



V / 



+ 



x 



\ 



2xy — y 2 




Bb 
Pi 



J 



\ 



y (y 2 - 3s 2 ) 
-x (x 2 — 3y 2 ) 




(40) 



/ 



These are the quadrupolar, hexapolar and octopolar components of the field. As long as 
Po ^> 1 one can neglect the higher order terms ending up with equation (JTJ) for B = 0. The 



gradient b = —^f - is fully determined by the strength Bb of the bias field and the current /. 

For our discussion we choose the experimental parameters I = 2 A and Bb = 1CT 2 T. 
These are typically achievable values which are used for realizing magnetic microtraps. For 
this setup one finds the line of vanishing field strength at a distance po = 40pm above the 



wire. The gradient b evaluates to b 
(2 2 Si,F = |) one obtains g = \gp\ 



250—. For the Lande-factor of 6 Li in its ground state 



| — ~|. Together with the mass M = 10964. 67m e the 
energy and the length scale become ^ = 0.864 • 10~ 9 eV and (^) _ ^ = 89.63nm, 

respectively. The resonance energy of the ground state is E = l.YlneV which corresponds 
to a temperature of 12.97 pK. This temperature regime is certainly accessible by todays cold 
atom experiments. Correspondingly the transition frequencies to excited states lie in the 
MHz regime. The lifetime of the ground state evaluates to 2.04/is. Under consideration of 
the exponential scaling with increasing angular momentum almost arbitrarily long lifetimes 
can be achieved by preparing the atoms in high A 2 eigenstate, e.g. the minimum lifetime 
for a state with |m| = y is 46.61s. 



VI. CONCLUSION AND OUTLOOK 

We have investigated the motion of a neutral spin-|-fermions in a magnetic quadrupole 
guide. The impact of an additionally applied homogeneous Ioffe field is also studied. Intro- 
ducing a canonical scaling transformation of the phase space coordinates we have derived 
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an effective two-dimensional Hamiltonian depending on a single parameter 7. The energies 
and decay widths of resonance states of the Schrodinger equation have been calculated by 
employing the complex scaling method. Utilizing a 2-dimensional harmonic oscillator basis 
we were able to converge hundreds of resonance states. 

The analysis of the underlying Hamiltonian revealed a large number of symmetries. In 
the absence of a homogeneous Ioffe field we have found 15 discrete symmetries of both 
unitary and anti- unitary character in addition to the conserved quantity A z = L z — S z . A 
deeper investigation of the underlying symmetry group revealed a two-fold degeneracy of 
any energy level. If a Ioffe field is applied A z is still conserved but only 7 discrete symmetries 
remain. Due to the altered symmetry group the degeneracies are lifted. 

We have calculated the resonance energies for several values of the parameter 7. A 
comparison to the values obtained by Hinds and Eberlein (Ref. [^]) has been performed. For 
7 = (vanishing Ioffe field) the resonance energies and decay widths form a regular pattern 
in the E — V plane. Increasing the Ioffe field leads to a distorted distribution and the decay 
widths are pushed towards larger values. Thus the stability of the resonance states increases 
with the Ioffe field strength. An analysis of the lowest resonance has shown an exponential 
increase of lifetimes which agrees with results obtained by others [l5j. Furthermore we 
could show an exponentially increasing lifetime with increasing modulus of the m quantum 
number. Apparently, with increasing \m\ the wavefunctions become localized farther from 
the center of the guide where transitions to continuum states are induced. 

We could show the existence of so called quasi-bound states which can be described by 
a scalar radial Schrodinger equation. The approximate eigenenergies agree very well with 
the resonance energies obtained from the complex scaling calculation and become exact in 
the limit of high m quantum numbers. But even for low angular momenta an astonishing 
agreement could be observed. For 7 = (without Ioffe field) this is due to the fact that 
m is half-integer valued which leads to a non-vanishing angular momentum barrier. This 
prevents the particle from entering the center of the guide. However, the coupling matrix 
element to the unbound states does only assume significant values for p — > i.e. the 
corresponding transitions are strongly inhibited. For 7 = we have also calculated an 
analytic expression for the ground state energy in each m subspace. For 7 7^ the quasi- 
bound energies are compared to those obtained from the so-called adiabatic approximations. 
We have shown that our approach is in general more accurate and reproduces in particular 
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underlying degeneracies. 

The results have been applied to the case of 6 Li in the F — | state. Here we have con- 
sidered a magnetic guide generated by a current carrying wire together with a homogeneous 
bias field. We have shown that for typical experimental parameter values the ground state 
energy corresponds to a temperature of a few micro-Kelvin. The lifetime of the resonance 
states can be extended up to minutes if the atoms are prepared in a sufficient high angular 
momentum state. 
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